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Abstract. It is shown that the endpoint of the first order transition line which merges into a crossover regime in the 
phase diagram of the Nambu-Jona-Lasinio model, extended to include the six-quark 't Hooft and eight-quark interaction 
Lagrangians, is pushed towards vanishing chemical potential and higher temperatures with increasing strength of the OZI- 
violating eight-quark interactions. We clarify a connection between the location of the endpoint in the phase diagram and 
the mechanism of chiral symmetry breaking at the quark level. Constraints on the coupling strengths based on groundstate 
stability and physical considerations are explained. 
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INTRODUCTION 



The determination of the phase diagram of QCD represents an interesting challenge both from the experimental and 
the theoretical points of view. On the experimental side several new facilities are aiming at an unprecedented look 
at the properties of hot and/or dense strongly interacting matter and from the theoretical point of view, as standard 
perturbative techniques fail in the low energy regime of QCD, alternative approaches have to be used. Great progress 
has been achieved over the last two decades in the theoretical study of the QCD phase diagram both in terms of 
effective low energy theories and lattice calculations (see for instance reviews [1,2, 3, 4] and paper [5]). 

The Nambu-Jona-Lasinio model [6, 7, 8], is regarded as an useful tool for the study of low energy hadron 
phenomenology as it shares with QCD the global symmetries and incorporates by construction a mechanism for 
dynamical chiral symmetry breaking (DjSB). The extension to include a 2Nf 't Hooft determinantal interaction 
(NJLH) [9, 10, 11, 12, 13] (where Nf is the number of flavors) breaks the unwanted Ua{1) symmetry exhibited 
by this model. The study of this model in the Ught quark sector (u, d and s) has shown a fimdamental flaw [14] (the 
absence of a globally stable ground state) which can be remedied with the inclusion of general scalar non derivative 
Ul{3) X Ur{3) symmetric eight quark interactions. 



FORMALISIVI 

The model Lagrangian that we are going to use, ^eff, can be decomposed in the following manner (see (1)): 
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The free Dirac term, J^d, contains a current mass matrix, m, diagonal in flavor space which explicitly breaks chiral 
symmetry (we wiU consider Mu = rrid ^ rris as such SU{3) flavor symmetry is also broken). The 4-quark NJL 
interaction, ^njl, if strong enough induces spontaneous chiral symmetry breaking (here Aq = -\/2/3 and the A^, 



a = 1, ... ,8 are Gell-Mann matrices acting in flavor space). The 6-quark 't Hooft term, ^h, breaks axial symmetry 
(here the Pl^r are chiral projectors and the determinant is over flavor space). The 8-quark interaction, can be spUt 

in two parts, one of which, -S?g^^ along with the 't Hooft term, induces OZl rule violation. 

This additional term ensures ground state stability as long as: > , +3^2 > and g\ > 1/G(k:/16)^ [14]. 
Through a suitable choice of parameters (a decrease in G with increased gi with remaining parameters kept fixed) 
the light scalar and pseudoscalar meson spectra in the vacuum can be left relatively unchanged, apart from a marked 
decrease in the a meson mass, and close to the one obtained in the NJLH model [15]. ^ 

This wide spectrum of allowed (G) encompasses different scenarios ranging from the low ^i(high G) case where, 
similarly to what happens in the NJLH model, the NJL term is the responsible for DjSB, to the high (low G) 
scenario the NJL term is too weak to induce it and the physical vacuum is in fact induced by the 6-quark interactions. 

These different scenarios can be disentagled by considering the non- vacuum properties of the model. This can be 
seen in study of the effects of the inclusion of a constant magnetic field [16, 17] or of finite temperature [18, 19] and/or 
chemical potential which are to a large extent dictated by the relative strength of the 4-quark and 8-quark interactions. 
Our focus will be on the latter since we are concerned with the phase diagram of QCD. 

For Nf > 2 ihe bosonization of this effective quark Lagrangian can be done in the path integral formaUsm with 
the aid of the introduction of two sets of flavour nonet bosonic fields: {(Ja,(l)a}, related to the physical scalar and 
pseudo-scalar mesons, and {sa,Pa} which are auxiliary fields [13]. The vacuum-to-vacuum amplitude can then be 
split in two parts, enabhng the separate evaluation of the functional integral over the fermionic degress of freedom and 
over the auxiliary variables. The former can be done using a symmetry preserving heat kernel scheme (generaUzed to 
deal with non-degenerate quark current masses [20, 21, 22]) whereas the latter can be done using the stationary phase 
approximation (for the present results it was done to leading order). 

The minimization of the effective (or thermodynamic) potential then corresponds to the self consistent solution of 
two sets of equations (for details see for instance [15] and references therein). One of these comes from demanding 
vanishing vacuum expectation value for the meson fields, the gap equations, 
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while the other comes from the Stationary Phase requirements to integrate out auxiliary fields {sa-,Pa}y 
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Here A, = M, — m;, with M;, the dynamical mass of the quark species i and the A,;, are differences of squared masses 
Aij = Mf — M?. The /, correspond to the averaged sum of /, = 1 /'iY,i=u.d.sJiy where the 7, are quark 1-loop euclidean 
momentimi integrals with i + 1 vertices. For their regularization we chose to use two Pauh- Villars subtractions resulting 
in the definition: 

|SA=l-(l-A'5^)/'5i' (4) 

Integrals of different orders can be easily related through recursion relations: = —dj^iJi- Their generalization to 
take into account the effect of finite temperature and/or chemical potential is done through the introduction of the 
Matsubara frequencies as per usual. 



' We chose to fit the model parameters G, K, g2. A, m„(= mj), to reproduce the following vacuum properties: M^i = 0.138GeV, M/j; = 
0.494GeV, M^, = 0.958GeV, M^^ = 0.980GeV, = 0.092GeV and = 0.1 17GeV. This sets the stability bounds on the value for the coupling 

strength 347GeV~* = Si ' < Si < sf ' = 20300GeV~*. One might naively expect another constraint coming from the requirement that M(j is a real 

positive but the critical value for that to be false lies outside the stability conditions: gf^ = 20933GeV~*. 



The themodynamic (or effective) potential can be obtained by integration of the gap equations, resulting in (see 
[23] for details): 

= ^ (8G {hi + hl + h])+ iKhuhdK + 3g,{hl+hl+h])+6{hl + h^,+h',)) 
+ ^{J-i{Mu,A,Hu,T)+J-i{Md,A,Hd,T)+J-i{Ms,A,Hs,T)). (5) 
The quark loop integral contribution, the J-i integrals are given by: 
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where the last subtraction corresponds to the integrals evaluated with M = 0. In the last line we have made explicit 
the separation between the vacuum and the medium contribution. The latter involves (anti-)quark occupation numbers 
which are defined as: 

- 1 _ 1 n^ 

RESULTS 

As was mentioned above there is a wide range of allowed values for ^1 (and G, as they have to be changed 
in consonance for the meson spectra to be left more or less unchanged) if we take into account only vacuum 
stabiUty restrictions. In [18, 19| it was shown that the temperature behaviour at vanishing chemical potential changes 
dramatically over this range. The OZl violating eight quark interactions dictate the temperature, slope and nature of 
the chiral transition: higher gi results in lower Tc, with a steeper slope and above a certain critical value one obtains 
a first order transition (which is made possible by the above described scenario of a physical vacuum induced by the 
't Hooft term in a potential that would otherwise correspond to a distorted Wigner-Weyl phase). As there is growing 
evidence coming from lattice calculations [24, 25] that the transition at = is crossover this can be interpreted as 
another upper bound for ^1 ^. 

In the case for zero temperature and finite chemical potential (considered equal for all quark species /i = /i„ = 
l^d = Ms)> on the other hand we always find a first order chiral transition for some critical value Hc which is lower for 
higher ^1. The analysis of the dependence of the quarks dynamical mass as a function of the chemical potential for 

(4) 

different values of gi (see Fig. 1) reveals that, below a certain critical value g\ , the model describes a low density 
phase of massive quarks (notice that on the left panel the line refering to M„^^(/i) crosses the diagonal before the first 
order transition jump). As this is deemed unphysical (see for instance [26]) this sets a lower bound for the interaction 
strength ^. 

Ploting the dynamical mass as a function on the bayonic density (as can be seen in the left panel of Fig. 2) we 
see furthermore that the jump in is higher and occurs at lower densities for stronger coupling. Furthermore the 
critical chemical potential for non- vanishing sttange quark density is lower for higher gi. Both these results may have 
interesting implications in the physics of compact stars. 

For finite T and jU (see Fig. 3) we find that the crossover/first order chiral transition line and the position of the 
Critical End Point (CEP) are highly sensitive to the chosen value for . As the transition as a function of T is steeper 
for higher gi the critical n for it to become fist order is lower for higher ^1: the CEP moves to the left in the phase 
diagram with increasing OZl-violating eight quark interactions. Furthermore in the studied cases the transition Une, 
Tc{li), for higher g\ lies completely inside that of lower gi. 



^ Using the already mentioned flttind procedure we obtain gf' = 8703GeV * which is lower than the upper bound required by vacuum stability. 
' For the chosen fits we have gf' = 2423GeV~* which is also inside the range allowed by vacuum stability 
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FIGURE 1. Dynamical mass of the quarks as a function of chemical potential at zero temperature {[M] — [fl] = GeV). Left panel 
refers to the weak coupling case set a) of Table 1 whereas the right panel refers to set b). The upper grey thick lines refer to Mj, 
the black thick lines refer to M„ and the dotted line refers to unstable solutions which are skipped in the process of the first order 
transition. The diagonal line refers to M = /x. 




a) b) 

FIGURE 2. a) Dynamical mass of the quarks {[M] = GeV) as a function of baryonic density (divided by po = 0. 17fm~^) at zero 
temperature. The upper grey and black lines correspond respectively to and M„^ in the weak coupling parameter set (the lower 
correspond to stronger coupling). The dotted lines correspond to unstable solutions which are skipped by the first order transition, 
b) Chemical potential above which the density of strange quarks becomes non- vanishing as a function of gi. The dashed line 
corresponds to the range of gi values alUowed by vacuum stability but excluded by the requirement of the unexistence of a low- 
density phase of massive quarks. The full line goes through the allowed values up to the limit that corresponds to the requirement 
of crossover transition at zero chemical potential. 

CONCLUSIONS 

The eight quark interactions necessary to stabiUze the ground state in three flavor NJLH model, dictate the phase 
diagram of the model. Their strength can be constrained using considerations on vacuum stability and the requirements 
of crossover at vanishing chemical potential as well as the imexistence of a low density phase of massive quarks. 



TABLE 1. Parameter sets obtained by fits to the ;;r , k, rj' and masses, as well as the n and K weak 
decay constants. The value for OZI violating part of the eight quark interactions (marked with *) is not fitted. 



Sets 


mu 
(MeV) 


Ms 

(MeV) 


(MeV) 


Ms 
(MeV) 


A 

(MeV) 


G 

(GeV-2) 


K 

(GeV-5) 


Si 

(GeV-^) 


g2 

(GeV-^) 




a 1 


5.9 


186 


359 


554 


851 


10.92 


-1001 


1000* 


-47 


b 1 


5.9 


186 


359 


554 


851 


7.03 


-1001 


8000* 


-47 



T 




FIGURE 3. Phase diagram for (liferent values of ^i, given in units of GeV^^ . Temperature, T, and chemical potential, /i, are 
given in GeV. The gray lines correspond the crossover transition for the two parameter sets from Table 1 (outermost corresponds to 
weaker coupling). The blue dots, the full and dotted black lines correspond (again for the two sample parameter sets) respectively 
to: CEP, first order transition and spinodals. The red dot marks the position of the CEP for the critical value below which a low 
density gas of massive quarks is obtained T = 0, gi = 2423GeV^**. The dashed line corresponds to the range of allowd CEP. 
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